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ABSTRACT 

Sirius has always attracted a lot of scientific interest, especially after the discovery of a 
companion white dwarf at the end of the 19*^ century. Very early on, the existence of 
a potential third body was put forward to explain some of the observed properties of 
the system. We present new coronagraphic observations obtained with VLT/SPHERE 
that explore, for the very first time, the innermost regions of the system down to 0.2" 
(0.5 AU) from Sirius A. Our observations cover the near-infrared from 0.95 to 2.3 /im 
and they offer the best on-sky contrast ever reached at these angular separations. 
After detailing the steps of our SPHERE/IRDIFS data analysis, we present a robust 
method to derive detection limits for multi-spectral data from high-contrast imagers 
and spectrographs. In terms of raw performance, we report contrasts of 14.3 mag at 
0.2", '^16.3 mag in the 0.4-1.0" range and down to 19 mag at 3.7". In physical units, 
our observations are sensitive to giant planets down to 11 Mjup at 0.5 AU, 6-7 Mjup 
in the 1-2 AU range and ^4 Mjup at 10 AU. Despite the exceptional sensitivity of our 
observations, we do not report the detection of additional companions around Sirius A. 
Using a Monte Carlo orbital analysis, we show that we can reject, with about 50% 
probability, the existence of an 8 Mjup planet orbiting at 1 AU. 

Key words: methods: data analysis - techniques: high angular resolution - star: 
individual: Sirius A - planetary systems 


1 INTRODUCTION 

As the brightest star in the sky, Sirius has attracted the 
attention of humanity since prehistoric times. After the dis¬ 
covery of a second component in the Sirius system in 1862 
(Bond 1862), 18 years after its theoretical prediction (Bessel 
1844), the system has been extensively monitored using var¬ 
ious observing techniques. 

In addition to the discovery of Sirius B as a white dwarf 
companion, there are particularities of the Sirius system 
that continue to generate much scientific interest. Sirius A 
and B form a close binary system, with a 50 year period 
(van den Bos 1960), on an elliptical orbit (e = 0.59) that 

* Based on observations made with ESO Telescopes at the La 
Silla Paranal Observatory under program IDs 60.A-9365 and 
60.A-9382. 
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brings Sirius B closer to Sirius A (8 au) than Saturn from 
the Sun. Sirius A is known to have the earliest spectral type 
amongst white dwarf companions (Holberg et al. 2013), and 
Sirius B has one of the highest known masses for a white 
dwarf (Gatewood & Gatewood 1978). 

Ancient Greek and Chinese records also suggest that 
the colour of Sirius may have changed within historical 
times (Gry & Bonnet-Bidaud 1990). One of the possible 
scenarios put forward to explain this change in colour is 
related to a transient reddening of Sirius caused by mat¬ 
ter ejected by tidal interaction with a potential third body 
(Bonnet-Bidaud & Gry 1991). Irregularities in the orbital 
motion of Sirius B were detected very early on suggesting 
the existence of a third component in the system. Detec¬ 
tion of possible visual companions (e.g. Baize 1931) and of 
a 6.3 years periodicity perturbation in the orbit of the Sir¬ 
ius A-B system (Benest & Duvent 1995) has fed hopes of 
detecting a low mass companion in the system. In particu- 
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lar, the periodic perturbation reported by Benest & Duvent 
(1995) led them to predict the existence of a companion of 
maximum mass ~50 Mjup orbiting Sirius A at less than 
3” (7.9 au). Because of the relatively young age of Sirius A 
(225-250 Myr; Liebert et al. 2005), a companion of this mass 
would be easily detectable with imaging in the near-infrared 
(e.g. Burrows et al. 1997; Baraffe et al. 2003). 

These particularities make the search for stellar or sub- 
stellar companions to Sirius A, or B, with imaging particu¬ 
larly attractive. However, the extreme brightness of Sirius A 
has long been known to represent a major issue for such 
a search. Several attempts at detecting possible compan¬ 
ions by imaging the close stellar field around Sirius from the 
ground using coronagraphic devices (Bonnet-Bidaud & Gry 
1991; Bonnet-Bidaud et al. 2000) and adaptive optics 
(ESO/ADONIS; Bonnet-Bidaud & Pantin 2008), or from 
space with the HST (Schroeder et al. 2000), have proven 
unsuccessful. Deep imaging performed with ADONIS in 
K band reached a limiting contrast of 9.5 mag at 3” to 
13.1 mag at 10” from Sirius A, this translates to the abil¬ 
ity to detect massive brown dwarf companions. The most 
recent and most sensitive study to date was performed by 
Thalmann et al. (2011) using Subaru/IRCS at 4.05 fim in 
saturated imaging, which allowed them to achieve a sensi¬ 
tivity to companions of 6-12 Mjup at 1”, 2-4 Mjup at 2” 
and 1.6 Mjup beyond 4”. While they provide very strong 
constraints against the existence of the massive companion 
with a 6.3-yr period predicted by Benest & Duvent (1995), 
their sensitivity to giant planets at very small angular sepa¬ 
ration (<0.7”, 1.8 au) is limited by the saturation of Sirius A 
on the detector. Thalmann et al. (2011) conclude their analy¬ 
sis with the hope that new-generation high-contrast imagers 
would provide tighter constraints at very small separations. 

Four years later, we report the first observations of Sir¬ 
ius A with SPHERE (Spectro-Polarimetric High-contrast 
Exoplanet REsearch; Beuzit et al. 2008), the new generation 
high-contrast imager for the Very Large Telescope (VLT). 
We use this new instrument to perform the closest and deep¬ 
est imaging of (sub-)stellar companions around Sirius A. 
Apart from the strong scientific interest of potentially de¬ 
tecting either a stellar or planetary-mass object around Sir¬ 
ius A, these observations play the Important role of demon¬ 
strating sphere’s capability to observe nearby regions to 
the brightest stars. 

In Section 2 we present in more details the instrument 
and the observations. Section 3 details the data reduction 
steps performed to accurately calibrate the data. In partic¬ 
ular, we present several useful procedures for the new com¬ 
munity of SPHERE users. In Section 4 we perform the high- 
contrast analysis of our data to try to reveal the presence of 
faint companions, and in Section 5 we detail a new procedure 
to derive reliable detection limits for SPHERE near-infrared 
data, comparing it to more a straightforward approach. Fi¬ 
nally, we discuss our results and their implications in Sec¬ 
tion 6, before detailing our conclusions in Section 7. 


2 OBSERVATIONS 

Sirius A has been observed as part of the SPHERE instru¬ 
ment science verification under ESO programme ID 60.A- 
9365 (P.I. Gry). The SPHERE, instrument, installed at the 


VLT (Beuzit et al. 2008), is highly specialized being ded¬ 
icated to the high-contrast imaging and spectroscopy of 
young giant exoplanets. It is based on the SAXO extreme 
adaptive optics (AO) system (Fusco et al. 2006; Petit et al. 
2014; Sauvage et al. 2014), which controls a 41 x 41 actuator 
deformable mirror, and 4 control loops (fast visible tip-tilt, 
high-orders, near-infrared differential tip-tilt and pupil sta¬ 
bilization). The common path optics include several stress 
polished toric mirrors (Hugot et al. 2012) to transport the 
beam to the coronagraphs and scientific instruments. Sev¬ 
eral types of coronagraphic devices for stellar diffraction sup¬ 
pression are provided, including apodized pupil Lyot coro¬ 
nagraphs (Soummer 2005) and achromatic four-quadrants 
phase masks (Boccaletti et al. 2008). The instrument is 
equipped with three scientihc subsystems: the differential 
imager and spectrograph (IRDIS; Dohlen et al. 2008), an 
integral field spectrograph (IFS; Glaudi et al. 2008) and the 
ZIMPOL rapid-switching, imaging polarimeter which works 
in the visible (Thalmann et al. 2008). 

Observations presented here were acquired on the 2014 
December 07 in the IRDIFS_EXT mode, where both the IFS 
and IRDIS observe in parallel. In this setup, the IFS covers 
wavelengths from 0.95 to 1.65 fim (Y-H bands) at a reso¬ 
lution of ~30 in a 1.7” x 1.7” square field of view (FoV), 
while IRDIS observes in the dual-band imaging mode (DBI; 
Vigan et al. 2010) with K12, a pair of filters in the K band 
(Aki = 2.110 nm, Ak 2 = 2.251 /rm, ~0.1 fim bandwidth), 
within a ~4.0” radius circular FoV. The observations were 
performed in pupil-stabilized mode to perform angular dif¬ 
ferential imaging (ADI; Marois et al. 2006) with an apodized 
pupil Lyot coronagraph designed for the K band (ALC_Ks), 
which uses a coronagraphic mask of radius 120 mas and an 
optimized Lyot stop. The IRDIS detector was dithered on 
a 4x4 pattern to reduce the effect of the residual flat-held 
noise. 

Originally the programme was allocated two hours to 
observe Sirius A and B, each for one hour. Unfortunately, 
the observations of Sirius B could not be performed due 
to the brightness of Sirius A, located 9.9”away at the time 
of the observations. During the observations, bright stray 
light diffracted by the telescope spider vanes entered the 
FoV of the wavefront sensor, causing an over-illumination 
and an AO loop crash, and prevented further observations 
of Sirius B. To take advantage of the full time allocation of 
the programme, the sequence on Sirius A was then extended 
after the failed attempt on Sirius B. 

The complete sequence on Sirius A spanned 2 h and 30 
min, totalling 106 deg of FoV rotation. The time was divided 
into four continuous sequences, with interruptions of 10, 23 
and 8 min in-between, the longest interruption correspond¬ 
ing to the attempt at observing Sirius B. For IRDIS, the 
total integration time was 88 min, divided into data cubes 
of DITxNDIT^=4x25 s, while for the IFS it was 59 min, 
divided into data cubes of DITxNDIT=2x34 s. The differ¬ 
ence in total integration time between the two instruments 
is only the result of detector readout overheads that are 
more important for the IFS (efficiency of 51% for DIT=2 s) 
than for IRDIS (efficiency of 77% for DIT=4 s). The con¬ 
ditions were excellent, with seeing values below 0.5”for the 

^ DIT: detector integration time; NDIT: number of DIT. 
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Table 1. Observing log on 2014 December 07 


Seq. 

LST time 
(hr:min) 

IFS 

DITxNDITxNEXP 

Texp 

(min) 

IRDIS 

DITxNDITxNdither 

Texp 

(min) 

FoV rot. 
(deg) 

Seeing®' 

(arcsec) 

Sr'^ 

(%) 

1 

04:41-05:13 

2x34x16 

18 

4x25x16 

27 

4 

0.52 ±0.10 

92.8 ±1.6 

2 

05:23-05:32 

2x34x4 

5 

4x25x4 

7 

2 

0.46 ±0.05 

93.9 ±0.5 

3 

05:55-06:29 

2x34x16 

18 

4x25x16 

27 

29 

0.35 ±0.05 

95.0 ±0.6 

4 

06:37-07:11 

2x34x16 

18 

4x25x16 

27 

50 

0.42 ±0.09 

94.5 ± 1.0 

all 



59 


119 

106 

0.43 ±0.10 

94.1 ± 1.4 


The seeing and Strehl ratio estimations are calculated over periods of 10 s every 30 s by the real-time computer. The Strehl ratio is 
expressed at 1.6 pm. The error bar is calculated as the standard deviation of the values. 


duration of the observations. The observing conditions as 
reported by the observatory were clear. Combined with the 
extreme brightness of Sirius A {R = —1.46), this resulted 
in exquisite AO performance with reported Strehl ratios of 
over 90% for all observations. For sequences 1, 2 and 4, the 
AO spatial filter was set to a size of 1.31 A/D, and for se¬ 
quence 3, during the best observing conditions, the size of 
the spatial filter was decreased to an even smaller size of 
1.14 A/D. The coronagraphic observations of Sirius A are 
summarized in Table 1. 

The time between long coronagraphic observations was 
used to acquire calibration data that will be used for the 
data analysis. These data include the following: 

• Stellar point-spread functions (PSFs) taken off-axis 
(~0.4") with the neutral density (ND) filter ND3.5, which 
provides the highest attenuation factors in SPHERE -- ap¬ 
proximately 3 300 in Y band, 16 300 in J band, 1200 in 
H band, and 700 in K band. For this calibration, the star is 
moved away from the coronagraph by applying an offset on 
the near-infrared differential tip-tilt plate. During this obser¬ 
vation, the AO visible tip-tilt and high-order loops remain 
closed to provide a diffraction-limited PSF. Three off-axis 
PSFs were acquired, one at the beginning of sequence 1 and 
two at the beginning of sequence 3. 

• ‘Star centre’ coronagraphic images where four symmet¬ 
ric satellite spots are created by introducing a periodic mod¬ 
ulation on the deformable mirror. This data is used in subse¬ 
quent analysis to determine an accurate position of the star 
centre behind the coronagraph, and hence the centre of field 
rotation. Four star centre frames were acquired, one during 
each of the observing sequences. 

• Sky backgrounds taken at a distance of ~30” from Sir¬ 
ius in a region of the sky empty of stars. These backgrounds 
are necessary for the IRDIS observations in K band, where 
the sky has a signihcant contribution. 

Data from the SPARTA real-time computer of the 
SAXO extreme AO system were collected at regular inter¬ 
vals in parallel of all the observations, including images from 
the differential tip-tilt sensor (DTTS). The DTTS removes 
a minute fraction of the incoming flux in the near-infrared 
(at 1.5 /rm) to image the PSF just before the coronagraph. 
This it used for the DTTS loop which keeps the PSF locked 
on the coronagraph once the observing sequence has started. 
These files also contain estimations of the seeing, wind speed 
and Strehl ratio made by the real-time computer. 

Standard calibrations for the IRDIFS_EXT mode were 
acquired in the morning as part of the instrument calibra¬ 


tion plan. For the IFS, this includes deep instrumental back¬ 
grounds with the same DIT values as the science data, de¬ 
tector flats at four different wavelengths and with a white 
lamp, an integral field unit (IFU) flat that allows to cali¬ 
brate the response of each individual lenslet, and finally a 
wavelength calibration frame where the IFU is illuminated 
using four different laser lines. For IRDIS, this includes deep 
instrumental backgrounds and a detector flat-held acquired 
in the K12 filter pair. 


3 DATA REDUCTION 

In this section we describe the data reduction of our Sir¬ 
ius observations, i.e. the steps that we followed to go from 
the raw data to cleaned and aligned data cubes that will 
later be used for our scientific analysis. Due to several issues 
with this data set, we choose to detail the data reduction 
quite extensively. We provide useful information that could 
be used for the reduction of SPHERE data by other users. In 
Appendix A, we describe some of the more technical details 
and we provide a link to our SPHERE/IFS pre-processing 
pipeline. 

We call pre-processing the steps of creating the basic 
calibrations and applying them to the science data. We dis¬ 
tinguish the case of IFS and IRDIS data because they are 
two conceptually very different instruments. For both instru¬ 
ments, we made use of the preliminary release (vO. 14.0-2) of 
the SPHERE data reduction and handling (DRH) software 
(Pavlov et al. 2008), as well as additional tools created for 
the purpose of improving or replacing the DRH. 


3.1 IFS pre-processing 

For the IFS, we followed the steps described in Mesa et al. 
(2015) for the creation of some of the basic calibrations using 
the DRH: backgrounds, master flat-field, IFS spectra posi¬ 
tions, initial wavelength calibration and IFU flat-field. After 
these calibrations, one would normally use directly the DRH 
to produce the (x, y, A) data cubes, but we introduce an ad¬ 
ditional pre-processing that performs the following actions: 

• accurate determination of the time and parallactic angle 
of each science frame based on header information (DATE-OBS 
and DATE FITS keywords for the start and end of each data 
cube); 

• subtraction of the background to each science frame; 

• normalization of the data by the value of the DIT and 
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by the attenuation of the ND filter used using their trans¬ 
mission curves^; 

• temporal binning of the raw science data to reduce the 
total number of frames. Sirius A, being extremely bright, 
required the minimum DIT of 2 s to be used, resulting in 
more than 1750 independent science frames. The binning 
was performed by directly averaging groups of four consec¬ 
utive frames. Given the excellent observing conditions and 
the overall quality of the data, we did not have to remove 
any bad frames, such as AO open-loops. The binning was 
performed within each of the data cubes, i.e. we did not bin 
frames coming from two different cubes. This resulted in 
the loss of ~150 frames over the whole sequence because the 
NDIT of the cubes was not a multiple of 4. After the bin¬ 
ning, the number of frames was reduced from 1750 frames to 
400. The maximum field rotation per binned frame is equal 
to 0.38° at transit. It is important to perform the binning of 
the science data at the level of the raw science data rather 
than at a later stage after multiple interpolations, so as to 
avoid the addition of interpolation noise. During the binning 
process, the timing and position angle of each binned frame 
is determined using the information of each individual raw 
frame; 

• correction of bad pixels identified using the master dark 
and master flat-field DRH products - this step is introduced 
as a substitute to the bad pixel correction provided by the 
DRH. This cosmetic step is particularly important due to 
the large number of isolated bad pixels in the IFS detector. 
We compared several methods and found that the MASKIN- 
TERP IDL procedure®, performing a bicubic pixel interpola¬ 
tion using neighbouring good pixels as reference, gave the 
best results; 

• cross-talk correction. The light passing through the IFS 
lenslets forms one PSF for each lenslet. However, part of 
the light at each PSF position is actually contaminated 
by spurious signal from adjacent lenslets. Moreover, a sec¬ 
ond effect for the IFS spectra is that spurious light coming 
from other wavelengths can also contaminate the light at 
a considered wavelength. We refer the reader to the work 
of Antichi et al. (2009) for a theoretical treatment of the 
SPHERE IFS crosstalk. The consequence of this cross-talk 
effect on the IFS calibrated data cubes is the presence of 
spurious light from different wavelengths in each spectral 
frame. This effect is particularly evident when looking at 
images with the satellite spots used to centre each frame. 
In this case, a double satellite spot can be seen especially 
when looking to wavelengths where the total flux from the 
star is very low. We have developed a two-step procedure, 
that is applied to the raw data, to solve this problem. The 
first step corrects the small scale cross-talk using a kernel of 
41 X 41 pixels that contains a Moffat function of the type 
l/(l-|-rjigt/b|j;.), where raist is the distance from the centre of 
the array (expressed in pixels) and 6fac is an appropriate nu¬ 
merical value that has been determined iteratively by trying 
to minimize the brightness of the secondary satellite spots 
present at particular wavelengths as described above. The 
data is first convolved with the kernel, and the convolved 
image is then subtracted from the original data. To avoid 


^ https://www.eso.org/sci/facilities/paranal/ii 
® http://physics.ucf.edu/“jh/ast/software.html 


subtracting the value of the central pixels, the 3x3 pixel 
central region of the array is set to zero. The second step con¬ 
sists of removing the large-scale cross-talk that corresponds 
to the instrumental background in 64 x 64 pixel boxes. This 
step is not necessary when subtracting a proper instrumental 
background as we do in the second step of the intermediate 
pre-processing. 

The partially calibrated, temporally binned, frames are 
then injected into the DRH recipe that creates the data 
cubes by interpolating the data spectrally and spatially. In 
this step, no background subtraction or bad pixel correction 
are performed, since these calibrations have already been 
made. 

As mentioned previously, the IFS detector is strongly 
affected by bad pixels. Even with a thorough initial cleaning 
of the raw science frames, the produced {x, y, A) data cubes 
still show remaining bad pixels, mostly in the low signal-to- 
noise ratio (SNR) wavelengths corresponding to atmospheric 
absorption bands. To finish the pre-processing of the sci¬ 
ence data, these bad pixels are first identified using a sigma¬ 
clipping routine, and then corrected using the MASKINTERP 
procedure. 

An important additional step that we have implemented 
is a correction of the wavelength calibration, as performed by 
the DRH, which reduces calibration errors on from ~20 nm 
(at worst) to 1.4 nm. As this procedure is useful for the ex¬ 
ploitation of the instrument by the community, we describe 
it in detail in Appendix A2. As a byproduct of the wave¬ 
length calibration step we also obtain the spatial rescaling 
factor that will later be used to analyse the data (see Sec¬ 
tion 4). 

3.2 IRDIS pre-processing 

Pre-processing of IRDIS data is much more straightforward 
owing to its simple optical concept. This part of the analy¬ 
sis is handled by the LAM-ADI pipeline (Vigan et al. 2012, 
Vigan et al. submitted), which has been updated to handle 
multi-wavelength data from IRDIS and the IFS. Each of the 
images in the coronagraphic observing sequences are back¬ 
ground subtracted and divided by the flat-field. Bad pixels 
are corrected using the bad pixel maps created by the DRH 
which replaces them with the median of neighbouring good 
pixels. Since the IRDIS DBI mode uses well-characterized fil¬ 
ters, the scaling factor between the two wavelengths is taken 
as the ratio of the central wavelength of the two filters. We 
also apply some temporal binning of the data by averaging 
groups of four consecutive frames, resulting in a reduction 
of the number of frames from 1300 to 306. The maximum 
field rotation per binned frame at transit is 0.52°. 

3.3 Off-axis reference PSF 

Sirius A is the brightest star in the night sky, making its ob¬ 
servation difficult without a coronagraph. The ND3.5 filter 
used for our observations (see Section 2) was not sufficient to 
completely avoid saturation of the off-axis PSF of Sirius A. 
The PSF was slightly saturated or out of the linear range 
in the first five and last four spectral channels of the IFS, 
ir#ti4sy/@f311St@fltdi!a PSF in the K12 filters was saturated in 
IRDIS over a diameter of ^\/D. 
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To solve this problem, we make use of an observa¬ 
tion of /3 Pictoris performed the day after, also part of 
the SPHERE science verification under ESO programme 
ID 60.A-9382 (P.I. Lagrange). This target was observed in 
the same instrumental configuration (IRDIFS_EXT) and in 
extremely similar conditions (seeing of 0.30” ±0.03”, Strehl 
ratio of 94.6% ± 0.5%, clear observing conditions). The off- 
axis PSF was acquired using the same ND3.5, bnt without 
any saturation (maximum number of counts at ~12 000 , well 
below the linearity limit of 20 000). The data for the j3 Pic¬ 
toris off-axis PSFs from the IFS and IRDIS is calibrated 
and pre-processed in exactly the same manner as described 
in Sections 3.1 and 3.2. Then, this proxy PSF is scaled in flux 
to properly represent the Sirius PSF. Our scaling procedure 
is described in Appendix A3. We estimate an uncertainty on 
the final scaled PSF of 0.06 mag for IRDIS, and -^0.2 mag 
for the IFS modified spectral channels. 

3.4 Centring of the coronagraphic data 

During the observations of Sirius A, two of the atmospheric 
dispersion correctors (ADCs) were not tracking, one in the 
visible and one in the near-infrared. There are two major 
consequences to this malfunction: ( 1 ) there is residual dis¬ 
persion at the level of the coronagraph, which means that the 
star is not well centred behind the mask at all wavelengths, 
and (2) this centring is evolving with time as the other ADC 
of each pair continues tracking while its counterpart does not 
move. The first effect is mitigated by the fact that SPHERE 
uses a DTTS in its near-infrared arm to compensate the 
differential refraction between the visible (where the wave- 
front sensing is performed) and the near-infrared. Due to 
the DTTS, whatever happens with the ADCs, the star will 
remain centred on the coronagraph at the wavelength of the 
DTTS. This is very important since it provides us with a 
reference wavelength for the centring. 

For the IFS, we verify that this is indeed the case by 
measuring the position of the centre at every wavelength 
for each of the three star centre frames that were acquired"^ 
on Sirius A. As expected, at around 1.5 /rm the centres are 
within 0.5 pixel (3.5 mas) of each other, while at 1.0 /rm 
they are separated by more than 2 pixels. In addition, a 
visual assessment of the coronagraphic centring shows it to 
be accurate at longer wavelengths. This is expected as the 
original centring is performed in the H band. We adopt as 
the reference wavelength the spectral channel number 35 
(A 35 = 1.58 ^m), and we recentre all the other channels 
relative to it, as described in Appendix A4. 

The analysis is much simpler for IRDIS, due to the dif¬ 
ferential refraction attenuating with increasing wavelength. 
The three centres determined using the star centre frames 
are all within 0.16 pixel (2 mas), which means that the star 
centring on the coronagraph remained very stable through¬ 
out the full observation, despite the ADC tracking problem. 
A visual assessment of the data shows that the absolute cen¬ 
tring of the star behind the coronagraph appears satisfac¬ 
tory, with no apparent decentring. For the IRDIS sequence. 

One was acquired at the end of observing sequence 1, one 15 
min later at the beginning of sequence 2, and one 1 h 15 min later 
at the beginning of sequence 3. 


we adopt as the centre for the whole sequence the average 
of the three centres calibrated from the star centre frames. 
No additional recentring was required for IRDIS. 


4 SPECKLE NOISE ATTENUATION 

After performing the calibrations described in previous sec¬ 
tions, we are able to process our data using classical high- 
contrast imaging analysis techniques. The goal of these tech¬ 
niques is to attenuate or remove the stellar halo and speckle 
pattern that limit the ability to detect faint companions. 
For this purpose, we use both the angular diversity pro¬ 
vided by ADI (Marois et al. 2006) and the spectral diversity 
provided by spectral differential imaging (SDI; Racine et al. 
1999). Over the years, many data analysis algorithms have 
been developed to take advantage of such diversity. We use 
as our primary means of investigation a principal compo¬ 
nent analysis (PCA) following the formalism described by 
Soummer et al. (2012). Because we have multi-dimensional 
data (temporal and spectral), there are several possibilities 
that can be used to build the PCA modes and subtract the 
speckle pattern. Below we describe the three ways we use to 
analyse our Sirius data set. 

4.1 Angular differential imaging analysis 

The first technique consists of looking at each spectral chan¬ 
nel of the IFS and IRDIS independently and only apply ADI 
to remove the speckle pattern. For any given spectral chan¬ 
nel, the PCA modes are constructed using all images of this 
channel, resulting in a total of 400 modes for the IFS and 306 
for IRDIS. The modes are calculated over the whole image, 
from an inner radius of 0.2”, up to 0.8” for the IFS and up 
to 3 . 7 ” for IRDIS. While not necessarily the absolute best 
approach in terms of performance, the quality and stability 
of the data ensures excellent results with this scheme, while 
keeping the computing time reasonable. An increasing num¬ 
ber of modes (up to 30% of the total number) is subtracted 
to each of the science frames, before applying derotation and 
averaging all the speckle-subtracted images. 

A variant of this method, that we also explore, is to bin 
the data spectrally before applying ADI. For the IFS we bin 
the data into the equivalent Y, J and H broad-band filters, 
and apply the same analysis as described above. For IRDIS, 
we simply add the K1 and K2 data to form a single image. 

These methods are equivalent to considering that we 
have independent data sets at different wavelengths. While 
it is efficient for the subtraction of the speckles, it does not 
take advantage of the spectral dimension to further reduce 
the level of speckle-noise. 

4.2 Angular and spectral differential imaging 

Our third method uses the spectral information to perform 
SDI in addition to ADI. The PCA modes are constructed 
using all the frames, after rescaling them spatially (see Ap¬ 
pendices A2 and B) in order to make the speckle pattern 
match at all wavelengths. In the rescaled images, the speckle 
pattern is stable while the signal of a companion is moving 
with time and wavelength. This is more aggressive in terms 
of speckle subtraction, but the large spectral range covered 
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Figure 1. Final images obtained after ADI or SDI+ADI analysis of the IFS and IRDIS data. Left: IRDIS FoV after ADI analysis in 
the K1 filter. Right, top row: IFS FoV after ADI analysis of the data spectrally binned to equivalent J and H broad-band filters. Right, 
bottom row: IFS and IRDIS inner FoV after the full SDI-I-ADI analysis. The speckles and stellar halo have been subtracted using a PCA 
analysis (see text for details about the analysis). 


by the IFS (Y-H) gives room for using SDI with a safe mar¬ 
gin. The bifurcation radius (Thatte et al. 2007) is equal to 
148 mas. This defines the separation above which the signal 
of a companion will have moved by more than A/D after the 
spatial rescaling. Above this radius we will be able to build 
a reference PSF that is not significantly biased by that of 
the companion. With this approach, we obtain 15 600 modes 
for the IFS data set, and 612 modes for IRDIS. Here again, 
the modes are calculated over the whole image. For the IFS 
analysis, we subtract up to 5 000 modes in steps of 50 modes, 
and for IRDIS we subtract up to 500 modes in steps of 10 
modes. After the modes subtraction, all frames are spatially 
rescaled back to their original size, derotated and averaged 
to obtain the final image. 


4.3 Results 

The three analyses are performed from an inner radius of 
0.2" for both the IFS and IRDIS. This inner radius is not 
extended down to the theoretical inner-working angle of the 
coronagraph, 120 mas, due to the ADC tracking problem. 
After the recentring of the frames with respect to the real 
star centre (see Section 3.4), the coronagraphic mask is offset 
from the centre of the images. We estimate that, in the worst 
case, the coronagraphic mask covers up to 0.18”. We add a 
small margin and provide quantitative measurements from 
0 .2" outwards. 

All the final images were inspected visually and com¬ 
pared to each other to look for point-like structures. The 
multi-wavelength coverage was used to discriminate any 
speckle residual from the signal of any potential real detec¬ 
tion. We have not identified any point-like structure in any of 


the images and all visible structures in the images have been 
identified as residual speckle noise by our multi-wavelength 
comparison. Fig. 1 presents several of the final images ob¬ 
tained with either the IFS or IRDIS. Although some struc¬ 
tures are still visible after subtraction of the speckles, none 
of them can be identified as a real detection. 


5 DETERMINATION OF THE DETECTION 
LIMITS 

5.1 Method description 

Deriving detection limits in high-contrast imaging is getting 
increasingly difficult as the new generation of instruments 
incorporate high-performance dual-band imagers, such as 
IRDIS, or IFSs. The determination of detection limits in 
DBI data in the case of no-detection has already been stud¬ 
ied previously in the context of the SPHERE/IRDIS scien¬ 
tific preparation or surveys with VLT/NACO in SDI mode 
(Vigan et al. 2010; Maire et al. 2014; Rameau et al. submit¬ 
ted). There are essentially two approaches that have been 
explored so far: (1) the semi-analytical approach where the 
residual flux after SDI reduction is predicted as a func¬ 
tion of separation and physical parameters of a planet using 
the evolutionary model tables, allowing a direct conversion 
of the residual noise in the speckle-subtracted images into 
mass/temperature (Rameau et al. submitted); (2) the data 
analysis approach where one introduces in each image fake 
planets with flux based on simulated spectra (Maire et al. 
2014). Whatever the approach, it is generally agreed that 
the result of data analysis involving SDI cannot be directly 
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converted to a physical mass limit because of the degeneracy 
induced by the spectral difference. 

For IFS data, the relative novelty of this type of instru¬ 
mentation in high-contrast imaging has not yet triggered 
extensive studies about how to best derive accurate physical 
detection limits. This is however a very important problem 
as all large direct imaging surveys that recently started us¬ 
ing IFSs will require the derivation of such limits for a large 
number of stars to be able to perform the hnal statistical 
analysis of the giant planet population. 

We present here a method that we have developed for 
VLT/SPHERE IRDIFS data, based on the second approach 
described above for SDI data, but modified for our simul¬ 
taneous IRDIS and IFS data. In summary, we inject fake 
planets into the data, perform the speckle removal analysis, 
and then measure the level at which each of the planets is 
detected to establish a detection limit. To make the analysis 
more efficient, several fake planets are injected simultane¬ 
ously in each analysis at the same contrast level. They are 
placed at increasing position angles and separation from the 
star to form a spiral. Increments in position angle and sepa¬ 
ration are chosen so as to make sure that the self-subtraction 
effects, created radially by SDI and azimuthally by ADI, do 
not overlap for different fake planets. The fake planet pat¬ 
tern is identical between the IFS and IRDIS, and for each 
contrast level the pattern is introduced at 10 different orien¬ 
tations to improve the final detection statistics. The smear¬ 
ing induced by the FoV rotation in each of the frames is 
properly taken into account when creating the fake planets 
map. 

For the speckles subtraction, we perform for each newly- 
created data set the same analysis as described in Section 4 
using the SDI-(-ADI approach. At the end of each analysis, 
the final image is convolved with a A/D aperture, in order 
to remove noise at spatial scales smaller than the expected 
planet signal. Then we define the SNR of each of the injected 
planet as the maximum value of the convolved image at the 
known location of the planet divided by the rms of statisti¬ 
cally independent pixels (i.e. separated by more than A/D) 
in an annulus located at the same separation as the planet 
in the convolved image. Within this annulus, we exclude the 
signal of the planet itself, the area located at ±25 deg on 
either side that contains significant self-subtraction effects, 
and the radial spectral self-subtraction effects from other 
fake planets located at other separations and position an¬ 
gles. 

We perform two distinct analyses that serve different 
purposes. In the first, detailed in Section 5.2, we introduce 
fake planets with a stellar-like spectrum, i.e. with the same 
contrast ratio with respect to the star in each spectral chan¬ 
nel. In the second, detailed in Section 5.3, the contrast of 
the fake planets is based on synthetic spectra representing 
planets at different masses. 

Our analysis focuses only on the detection of compan¬ 
ions, as opposed to their characterization. Indeed, these are 
two distinct aspects of the same problem that do not require 
the same treatment. While the data analysis focused on char¬ 
acterization depends on the accurate calibration of all the 
biases on the photometry and astrometry, data analysis fo¬ 
cused on detection can afford to be much more aggressive, 
at the expense of quantitative flux measurements. 



Angular separation [as] 

Figure 2. SPHERE sensitivity in contrast to obtain a detec¬ 
tion at SNR=5 with the IFS (plain curve) and IRDIS (dashed 
curve) assuming a companion with a constant contrast ratio with 
respect to the star in all spectral channels. We show the sensi¬ 
tivity obtained both for an SDI-j-ADI analysis and for an ADI- 
only analysis (see text for details). The self-subtraction effects at 
small separation in IRDIS DBI mode strongly affect the perfor¬ 
mance below 0.7", while the wide spectral coverage of the IFS 
(0.96-1.6 pm) allows to mitigate this effect and explore smaller 
separations. These limits are compared to the limits published by 
Thalmann et al. (2011) with Subaru/IRCS in the Br a filter at 
4.05 pm. 

5.2 Raw contrast detection limits using fake 
planets with stellar-like spectra 

Introducing fake planets with a stellar-like spectrum, i.e. 
with a constant contrast ratio with respect to the star, is 
not physically realistic, but it is useful to get an idea of the 
raw performance, which can be compared between instru¬ 
ments, since it is in no way model-dependent. It is, however, 
pessimistic when using SDI because the self-subtraction ef¬ 
fects at small angular separation will significantly affect the 
sensitivity. 

Fig. 2 shows the result of this analysis for the IFS and 
IRDIS. The detection limit corresponds to the limit where 
fake planets are detected at an SNR of 5 or higher. The 
penalty term induced by small sample statistics is taken into 
account (Mawet et al. 2014). The IFS detection limit shows 
excellent performance, with a contrast of 14.30 mag reached 
at 0.2" and 16.3 mag at 0.4". There is a slight increase of the 
limit around 0.6”, which corresponds to the signature of the 
AO-cutoff frequency averaged over all wavelengths, as all 
wavelengths are combined together for this analysis. This 
represents the best contrast ever reached in high-contrast 
imaging at these angular separations. 

For IRDIS, the fact that only two wavelengths are avail¬ 
able significantly affects the performance at small separation 
due to the spectral self-subtraction effects. There is a steep 
rise of the limit for separations <0.5”, but the limit joins 
that of the IFS at 0.6”, and it remains rather flat around 
16.7 mag in the 0.8-1.1” range. Similarly to the IFS, the AO- 
cutoff frequency in K band creates a slight increase around 
1.1”. Then the increase in sensitivity follows nicely the shape 
of the halo to reach a contrast of 19 mag at 3.7”. In Fig. 2 we 
compare our SPHERE limits with the best limits on Sirius A 
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published by Thalmann et al. (2011), which were obtained 
in saturated imaging with Subaru/IRCS in the Br a narrow- 
band filter (4.05 /rm). In the speckle-limited regime of these 
observations (<3"), the gain of extreme AO and an efficient 
coronagraph becomes evident, with a gain of 3.5 mag at 1”, 

4.3 mag at 0.7” and access to separations as small as 0.2”. 
We would likely have accessed even smaller regions if the 
ADCs had worked properly. Nonetheless, these results show 
the huge potential of high-contrast imaging with the new 
generation of instruments. 

We also compare our SDI-I-ADI analysis with an ADI- 
only analysis for both the IFS and IRDIS. For the IFS, we 
bin the data spectrally to an H broad-band equivalent be¬ 
fore performing the ADI analysis, and for IRDIS we simi¬ 
larly add together the data of the K1 and K2 filters. The 
SDI+ADI performance for the IFS shows a clear gain of al¬ 
most 2 mag at 0.2” and ~1.5 mag beyond 0.4” with respect 
to an ADI-only analysis. There is no clear evidence that SDI 
is failing at these separations thanks to the very large band¬ 
pass of the IFS in the IRDIFS_EXT mode (52%), but the effect 
would surely become much stronger at smaller separations 
where we would get close to the bifurcation radius (148 mas, 
see Sect. 4.2). Another important point to keep in mind is 
that the use of PCA over temporal and spectral dimensions 
mitigates the adverse effect that SDI can have at small sep¬ 
arations, especially if the conditions are stable in time and 
the amount of field rotation is large. For IRDIS, we see that 
ADI starts to give a slightly improved sensitivity at 0.3” and 
below, but inside 0.3-2.0” the gain of using SDI is clearly 
visible. At larger separations the gain is only marginal. 

5.3 Physical detection limits using fake planets 
with realistic spectra 

The previous approach is useful to obtain an idea of the 
raw instrument performance in terms of contrast, but it 
is not realistic in terms of defining what physical objects 
can actually be detected. Indeed, the spectra of cool plan¬ 
etary mass objects are very far from flat, presenting strong 
molecular absorption features in the near-infrared when 
Tefj decreases (e.g. Allard et al. 1997; Burrows et al. 2006). 
Atmospheric model grids coupled with evolutionary mod¬ 
els (Burrows et al. 1997; Chabrier et al. 2000; Baraffe et al. 
2003) have been commonly used in imaging surveys or for in¬ 
dividual targets to convert single-band contrast limits into a 
mass limit a posteriori (e.g. Vigan et al. 2012; Meshkat et al. 
2013; Chauvin et al. 2015). However, multi-wavelength data 
obtained with an IFS or a dual-band imager adds a level of 
complexity. Spectral features have a strong impact on the de¬ 
tectability of these objects and the spectral self-subtraction 
cannot be calibrated easily when a large number of spectral 
channels are involved. 

For our analysis, we make use of the most recent the¬ 
oretical developments by using the BT-Settl atmospheric 
models (Allard et al. 2013; Baraffe et al. 2015), linked to the 
evolutionary models for the AMES-Cond grid (Baraffe et al. 
2003). The BT-Settl model couples a cloud model to the ID 
general stellar atmosphere code PHOENIX (Allard et al. 1994, 
2001; Hauschildt et al. 1997). The model considers the for¬ 
mation of a cloud deck composed of up to more than 200 
condensate species, with optical data for 50 species included 
in the radiative transfer. The grain size and density, and 



Figure 3. Comparison of the Sirius SED (Bohlin 2014) with 
BT-Settl atmospheric models (Allard et al. 2013; Baraffe et al. 
2015) calculated tor increasing masses from 2 to 20 Mjup using 
the Baraffe et al. (2003) evolutionary tracks. The dotted, dashed 
and dot-dashed vertical lines show the central wavelengths of 
the IFS spectral channels, IRDIS K1 and K2 DBI filters and Sub- 
raru/IRCS Br a filter respectively. The expected contrast for each 
planet is calculated per-channel by integrating numerically the 
flux of the planet model and dividing it by the flux of Sirius A. 


the abundances of chemicals in the gas phase, including the 
effect of element depletion induced by the grain formation, 
are computed layer per layer through the photosphere fol¬ 
lowing a comparison of the time-scales for nucleation, con¬ 
densation, gravitational settling or sedimentation, and mix¬ 
ing. Once rained out, the grains and their constituent ele¬ 
ments are removed from the composition of that layer and 
all higher levels of the atmosphere. These models can pre¬ 
dict the flux at the surface of a given object only defined by 
(log g, Teff), and the assumed (solar) composition. For the 
masses probed here, below ~15 Mjup the models predict 
the settling of most of the primary iron and silicate clouds 
below the photosphere. For lower Tgfr, corresponding to 5- 
10 Mjup, the formation of a secondary cloud layer within the 
photosphere, sitting above the primary cloud layer and con¬ 
sisting of MnS and Na 2 S, then of KCl, NaCl and some ZnS, 
is predicted, confirming the studies of Morley et al. (2012). 
In the 2014 release of the models used here, opacities for 
all the above revised alkali cloud species have been added, 
but a specific grid is calculated to cover a lower domain of 
Teff because of the age of our target. This extends into the 
regime of water ice and ammonium hydrosulphide forma¬ 
tion tentatively confirmed in the latest observed Y dwarfs 
(Beamfn et al. 2014). 

The spectra are simulated on a grid of Teu and log g, 
covering the range 200-1200 K and 3.5-4.5 dex respectively. 
We use evolutionary tracks calculated for the upper limit 
age of Sirius, 250 Myr (Liebert et al. 2005), to calculate the 
expected Tgff, log g and radius for planets of mass from 2 to 
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Figure 4. SPHERE sensitivity in mass to obtain a detection at 
SNR=5 with the IFS (plain curve) and IRDIS (dashed curve), 
based on the injection of fake planets with spectra based on 
the BT-Settl synthetic spectra and evolutionary tracks (see Sec¬ 
tion 5.3 for details). These limits are compared to the contrast 
limits of Thalmann et al. (2011) in the Br a filter at 4.05 fim 
that we converted into mass using the same set of models as for 
our SPHERE analysis. We discuss the comparison between the 
Subaru/IRCS and SPHERE/IRDIFS limits in Section 6 


15 Mjup by steps of 1 Mjup, and from 15 to 25 Mjup by steps 
of 2.5 Mjup. Then, for each wavelength from 0.9 to 2.3 fim, 
we interpolate into the (refr,log g) grid of synthetic spectra 
to produce the spectrum expected for each mass. Finally, to 
obtain the flux received on Earth (at the top of the atmo¬ 
sphere) from a planet orbiting Sirius, the spectrum is multi¬ 
plied by a dilution factor equal to {Rp/d)^, where Rp is the 
planet radius determined from the evolutionary tracks and 
d = 2.637 pc is the distance to Sirius. Then, the spectrum 
of each planet is numerically integrated in the IFS spectral 
channels and the two IRDIS filters to get an estimation of 
the integrated flux. To calculate the contrast at which the 
fake planets must be injected in each spectral channel, we use 
high-quality spectral energy distribution (SED) of Sirius A 
published by Bohlin (2014). The SED is also numerically 
integrated in each channel, and the contrast is simply ob¬ 
tained from the ratio of the planet flux with the Sirius flux. 
The Sirius SED and model spectra calculated for different 
masses are shown in Fig. 3. From this analysis, we expect 
contrasts in J band of 15.5 mag and 13.7 mag for 5 and 
10 Mjup planets respectively, and 17.5 mag and 14.3 mag 
for the same masses in K band. 

To estimate the detection limits, the process detailed 
in Sections 5.1 and 5.2 is repeated. Fig. 4 shows the final 
detection limits in physical units of projected orbital semi¬ 
major axis and mass of Jupiter. Here again the effect of small 
sample statistics is taken into account. Despite the extreme 
brightness of Sirius A, we are able to reach the planetary- 
mass regime with both the IFS and IRDIS inside of 0.5” 
(1.3 au). We are sensitive to masses below 10 Mjup for all 
projected semi-major axes outside of 0.7 au. Here again, 
the wide spectral coverage of the IFS shows its strength by 
allowing the detection of masses below 12 Mjup within the 
0.5-2.1 au range. In fact, the IFS detection limit is even 
better than IRDIS at the edge of its FoV. Similarly to the 
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Figure 5. Comparison of the IFS detection limits obtained us¬ 
ing our approach, based on the injection and redetection of fake 
planets with a realistic spectra, and using the conversion of the 
contrast limit from Fig. 2 using the evolutionary models in J band 
or in H band, as proposed in Mesa et al. (2015) and Zurlo et al. 
(submitted). 


stellar-like spectrum cased detailed Section 5.2, the IRDIS 
limit gets worse towards very small separations because of 
self-subtraction effects. 


5.4 Comparison to a simplified approach 

When the stellar SED is readily available or easy to de¬ 
termine by model-fitting (e.g. Vigan et al. submitted), the 
approach we have described in the previous sections is ex¬ 
tremely straightforward to implement. In particular, it has 
the advantage of solving one of the usual biases related to 
ADI and SDI analysis, which is the estimation of the self¬ 
subtraction effects that are not always easy to quantify, espe¬ 
cially when considering multi-wavelength data. Indeed, con¬ 
trary to the a posteriori conversion of detection limits usu¬ 
ally employed up to now, our approach embeds the model- 
dependent evolutionary part directly into the analysis of the 
data. At the end of the analysis, the detection level of the 
injected fake planets provides the mass detection limit di¬ 
rectly; therefore, the flux losses do not need to be deter¬ 
mined. The data need to be fully reprocessed many times 
to obtain a robust estimation of the detection limit which 
makes this process very computationally demanding. 

We use the Sirius data set to compare our approach with 
the less demanding one detailed in Mesa et al. (2015) and 
Zurlo et al. (submitted) for the estimation of the detection 
limit in physical unit. In their approach, the data analysis is 
also based on PCA, but the detection limit in contrast is es¬ 
timated from the residual noise in the final images, corrected 
for the self-subtraction effects using injection of fake plan¬ 
ets close to the noise floor, and normalized to the median 
peak flux of the PSF at all wavelengths. Then, this detec- 
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tion limit is converted into physical units using evolutionary 
models calculated in either the J- or H band broad-band 
filters. This approach provides faster results, as it is only 
necessary to fully process the data twice, however, it does 
not take into account all of the spectral information available 
in the IPS data. This informatioir is lost (at least partialy) 
when converting the final contrast detection limit using only 
predictions in a broad-band filter. 

The results of our comparison are showed in Fig. 5. In 
this figure we compare our IPS mass detection limit from 
Fig. 4 to the contrast limit from Fig. 2 that we converted 
into mass using the same set of models. To obtain the limit, 
the expected contrast between Sirius A and planets at in¬ 
creasing masses is calculated in broad-band J and H filters. 
Then, the final mass limit is obtained by interpolating this 
expected contrast at the values defined by the contrast limit. 
This approach shows very good agreement for the conversion 
using the H band, with differences of 5% at most with re¬ 
spect to the limit derived using the injection of fake planets. 
In J band the conversion is less accurate, with differences of 
10 % on average. 

The good agreement between our method and the con¬ 
version using the evolutionary models in the H band could 
be anticipated because the contrast is expected to be smaller 
in the H band: if a planet is detected in the IPS data, most 
of the strength of its signal will come from the H band; 
so converting the contrast detection limit using luminosity 
predictions in the H band will likely lead to a good approx¬ 
imation of the mass limit. Still, this comparison is an im¬ 
portant step for on-going large direct imaging surveys, as it 
shows that mass limits could be derived from simple contrast 
curves rather than going into much more advanced and time- 
consuming analyses. However, we note that one specific case 
does not make a rule. A more systematic study taking into 
account stars of different spectral types (and thus colours) 
and planetary atmosphere models would be required to draw 
more general conclusions. 


6 DISCUSSION 

In Section 5 we demonstrated that the new VLT/SPHERE 
instrument keeps its promises in terms of reachable contrast 
in the near-infrared for very bright stars. Fig. 2 illustrates 
the new parameter space, in angular resolution and contrast, 
that the new high-contrast imagers will access compared to 
the previous generation instruments. 

In the case of Sirius, the brightest star in the sky, this 
high-contrast at small separation translates to being sensi¬ 
tive to planets below 10 Mjup at separations of less than 1 au 
at its nominal age of 250 Myr. However, it is known from 
previous studies (Rameau et al. 2013; Skemer et al. 2014) 
that L band imaging is even more favourable for young gi¬ 
ant exoplanets detection because of the decreased stellar flux 
and increased planet flux at longer wavelengths. This is illus¬ 
trated in Fig. 3, where we see that the Sirius flux drops by a 
factor ~100 between J band and the Br a filter at 4.05 fim, 
and by a factor ~10 between K band and this filter. In this 
same range, the peak flux of low mass planets also increases 
by factors of 10 or more, resulting with expected contrasts of 
only 12.4 and 11.0 mag for 5 and 10 Mjup planets at 4.05 fim 
respectively. 



1 10 
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Figure 6. Mean probability of at least one detection of a sub- 
stellar companion around Sirius A in the combined Subaru/IRCS 
and SPHERE/IRDIFS data sets as a function of the companion 
mass and semi-major axis. We consider [a] the case of coplanar 
orbits with the Sirius A—B system {i = 136.6 deg) or [b] the case 
of unconstrained inclination i. Compared to the similar study by 
Thalmann et al. (2011), we are able to place constraints down to 
1 au for giant planets more massive than 7-8 Mjup thanks to 
our new SPHERE/IRDIFS observations, which offer sensitivity 
to planetary-mass objects down to 0.2” (0.5 au). 

For comparison of our IRDIFS detection limits to the 
Subaru/IRCS limits from Thalmann et al. (2011) in the 
Br a filter, we converted their contrast limits into mass using 
the same set of models. The use of this more recent grid of 
models has the effect of improving slightly their mass limits 
with respect to their published ones, which were obtained 
using the COND evolutionary models (Baraffe et al. 2003). 
The reader should however keep in mind that predictions in 
the 1-3 Mjup range are still very uncertain, whatever the 
type of model used. Still, using more advanced models is 
perfectly justified despite the uncertainties, which equally 
apply to the COND atmosphere models. This comparison 
of limits is particularly interesting, as it clearly shows the 
usefulness of having access to longer wavelengths for plan¬ 
ets detection. Indeed, Subaru/IRCS allows to reach lower 
masses than SPHERE: at 0.7” (1.8 au), they are sensitive 
to planets with masses as low as 6 Mjup (against '^7 Mjup 
with SPHERE), and they reach 1.5 Mjup at 1.5” (4 au). 

Despite these impressive limits in L band, one has to 
keep in mind the limitations of observing at long wave- 
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lengths: (i) the final sensitivity is usually dominated by the 
sky background noise except in a few very bright cases like 
Sirius A, and (2) the inner-working angle is also usually 
wider because of the larger wavelength, although new de¬ 
vices such as the annular groove phase mask (Mawet et al. 
2005) can provide access to inner-working angles almost 
down to the diffraction limit of 8 m telescopes (Absil et al. 
2013). Nonetheless, our comparison shows the exemplary 
complementarity between new high-contrast imagers in the 
near-infrared and other imaging instruments in the L- or 
M band (e.g. NACO at the VLT, NIRC2 on Keck). 

We use this complementarity to place constraints on the 
presence of additional companions around Sirius A. For this 
purpose, we perform orbital Monte Carlo simulations where 
planets are simulated on a grid of mass and semi-major axes, 
in the range 1-30 Mjup and 0.2-100 an respectively. At each 
point of the grid, 10® orbits are simulated with random ec¬ 
centricity, longitude of periastron, longitude of ascending 
node and time of passage at periastron. For the inclination, 
we simulate two cases: either a random inclination (statisti¬ 
cal weight proportional to sin i) or a fixed inclination equal 
to 136.6 deg, which corresponds to a co-planar case with the 
Sirius A-B system. For each simulated orbit, the expected 
projected separation of the planet is calculated for the three 
epochs where we have detection limits (two for IRCS, one 
for SPHERE), and we check the detectability of the planet 
with respect to the detection limits. For each point of the 
grid, the likelihood of detection is calculated as the fraction 
of planets detected at least one time over the total number 
of simulated planets. The results are presented in Fig. 6. 

Compared to the similar study by Thalmann et al. 
(2011), our new SPHERE/IRDIFS observations bring strong 
constraints on the existence of giant planets within the 0.6- 
2 au range where previous observations lacked good sensitiv¬ 
ity. The co-planar case is better constrained towards small 
semi-major axes, and in such a case we are able to reject 
with 90% probability the existence of a 10 Mjup planet at 
1 au. In the more general case this probability decreases 
to 70%, which is still high. At smaller separations or lower 
masses, the completeness values drop significantly and we 
cannot draw any firm conclusions. This is nonetheless the 
first time that we can probe with high confidence the re¬ 
gion below the critical semi-major Oc = 2.17 au where the 
simulations by Holman & Wiegert (1999) predict that sta¬ 
ble orbits can be found for planetary companions. Of course, 
we cannot reject the possibility of the existence of planets 
that could not be detected, either because their mass is too 
low to be detected by Imaging or because they are currently 
hidden behind or too close to the star. At small semi-major 
axes, the probability of the latter increases significantly and 
multiple epochs would be necessary to rule out giant planets 
at 1 au or below at high confidence. 

Finally, we conclude by noting that because Sirius A is a 
relatively massive young star of type AIV, it is not adapted 
to the search of low-mass companions by the radial velocity 
methods, so imaging remains one of the most viable options 
to detect a planet in this system. 


7 CONCLUSIONS 

The Sirius system is certainly among the most intriguing 
stellar systems at close proximity to Earth. We have used 
high-contrast imaging in the near-infrared with the new 
SPHERE instrument, recently commissioned at the VLT, 
to get the closest ever look at Sirius A. This new generation 
instrument combines extreme adaptive optics and modern 
coronagraphs to allow the detection of very faint planetary 
companions around nearby, young stars. 

As for any new instrument, especially when relatively 
complex like SPHERE, the data reduction of the first data 
sets by the community can be complicated. We have pro¬ 
vided some useful guidelines on how to improve basic cali¬ 
bration steps for the IFS, as well as how to deal with prob¬ 
lems that other users could encounter in the future such as 
saturated off-axis PSF and inaccurate or unstable centring 
of the star behind the coronagraph during the observations. 
Thanks to these improved calibration steps, we were able to 
get the most of this high-quality data set. 

In terms of raw contrast, we present what is the best 
contrast ever obtained in the 0.2”-0.5" range with a direct 
imaging instrument. Using fake planets with a stellar-like 
spectrum, we demonstrate that we reach with the IFS a 
contrast of 14.3 mag at 0.2" and 16.3 mag at 0.4". With 
IRDIS, the self subtraction effects limits our sensitivity at 
small separations (< 0.7"), but we demonstrate the very nice 
complementarity between the two SPHERE near-infrared 
instruments: IRDIS reaches the same sensitivity as the IFS 
at the edge of the outer IFS FoV, providing a perfect wider- 
field complement up to 4" where it reaches a contrast of 
~19.0 mag. This exceptional detection limit likely represents 
SPHERE at Its very best because of the very high brightness 
of the guide star [R = —1.46), the close to perfect observ¬ 
ing conditions (~0.4" seeing) and very large amount of field 
rotation (106deg). 

Despite their exceptional sensitivity, our observations 
bring another non-detection of companions around Sirius A. 
Multiple analyses of the IRDIFS data using approaches 
based on ADI or SDI-I-ADI do not reveal any reliable can¬ 
didate. Starting from this non-result, we presented a new 
approach for the computation of the IFS detection limits 
based on the use of the stellar SED, evolutionary models and 
associated atmospheric models. While recent high-contrast 
instruments are all equiped with an IFS, the problem of de¬ 
riving reliable detection limits for this type of high-contrast 
data has been mostly overlooked so far. However, it consti¬ 
tutes a necessary step for the large-scale direct imaging sur¬ 
veys that are currently on-going with SPHERE or its North- 
American equivalent GPL In this context, we compared our 
approach to a more simple one where the contrast detection 
limit is directly converted to mass using evolutionary pre¬ 
dictions in the H broad-band filter. The two approach are in 
good agreement in the Sirius case, however a more detailed 
study would be required to draw firm conclusions. 

Our limits around Sirius A show that we are sensitive to 
giant planets more massive than 11 Mjup at 0.5 au, 6-7 Mjup 
in the 1-2 au range, and down to ~4 Mjup at 10 au. Outside 
of this range, previous limits obtained by Thalmann et al. 
(2011) with Subaru/IRCS in the Br a filter at 4.05 /rm are 
more sensitive thanks to the decreased stellar flux and in¬ 
creased planet flux at longer wavelengths. This illustrates 
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the complementarity that exists between new near-infrared 
high-contrast imagers and L- or M band AO imagers for the 
search of planetary-mass companions. 

Monte Carlo simulations based on the 
SPHERE/IRDIFS and Subaru/IRCS detection limits 
show that our observations provide for the Erst time 
meaningful constraints on giant planets around Sirius A 
in the 0.5--1.5 an range, with about 50% probability for 
the existence of an 8 Mjup planet at 1 an. A second epoch 
observation with SPHERE would probably strengthen 
these conclusions because at small semi-major axes the 
probability increases that the planet is behind or too close 
to the star. In any case, additional SPHERE data will 
not bring better constraints for planets less massive than 
7-8 Mjup or at separations smaller than 0.5 au, so there is 
still plenty of room for smaller planets that could not be 
detected with direct imaging. The Sirius system will keep 
its cloud of mystery for the foreseeable future. 
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APPENDIX A: IMPROVING THE 
SPHERE/IRDIFS DATA REDUCTION 

A1 SPHERE/IFS pre-processing pipeline 

Hoping that it will help the community to work with 
SPHERE/IFS data, we make our pre-processing pipeline 
available at the following address: 

http://people.lam.fr/vigan.arthur/ 

It contains a set of tools that simplifies the process to go 
from the raw data to the properly calibrated (x, y, A) data 
cubes, providing a good starting point for subsequent SDI 
and ADI analyses. 


A2 Correction of the IFS wavelength calibration 

An important calibration for the IFS is the determination of 
the wavelength of each of the spectral channels in the data 
cubes. The DRH provides a good estimation of the wave¬ 
length of each channel, but one can easily realize that the 
accuracy is not sufficient by performing a spatial rescaling 
of each channel by a factor equal to Ao/Ai, where Aq and 
Xi are respectively the wavelength of the first and the ith 
spectral channels. After such a rescaling, all the diffraction 
features (AO control radius, halo, speckles) should remain 
stable when going through all the spectral channels. How¬ 
ever, when using the wavelength estimation from the DRH, 
the diffraction features are not completely stable after rescal¬ 
ing, which indicates that the wavelength estimation is im¬ 
perfect. To improve the wavelength estimation, we perform a 
new calibration in two steps that allow (1) to determine the 
relative scaling factor between the spectral channels, and 
(2) determine a reference channels for which the absolute 
wavelength is known precisely. 

For step (1), we use the star centre frames described 
in Section 2. In these frames, the periodic modulation intro¬ 
duced on the deformable mirror creates replicas of the stellar 
PSF that appear as four symmetrical satellite spots at a sep¬ 
aration of A/D from the star. The spacing of the spots 
is proportional to wavelength, so by measuring the spacing 
for channel i with respect to channel 0, we are able to deter¬ 
mine accurately the wavelength scaling factor. To increase 
the accuracy, we measure six different distances between the 
four spots (two between opposite spots, four as the sides of 
the squares). The final scaling factor is taken as the average 
scaling factor determined for the six distances, as illustrated 
in Fig. Al. There are small differences in the scaling factors 
determined independently, and from the departure between 
the average scaling factor and the independent factors, we 
estimate an accuracy better than 0.93% (0.26% on average) 
on the determination of the scaling factor. 

The scaling factor is only a relative measurement that 
does not provide information on the absolute value of the 
wavelength of each spectral channel. To improve the absolute 
wavelength calibration, we need to determine the wavelength 
of one reference spectral channel, and the wavelengths of all 
the others will then be extrapolated using the scaling factor. 


^ ^ Apjjroximate wavelength [pm] ^ ^ 

1.8n -1-'-^^^^^^^^^-r 



0 10 20 30 40 

Spectral channel index 


Figure Al. Scaling factor with respect to the first spectral chan¬ 
nel as a function of the spectral channel index. The top axis gives 
the approximate wavelength of the channels. In each channel, six 
distances are measured using the satellite spots to determine the 
average scaling factor (black curve). The inset image shows the 
star centre frame at a wavelength of 1.24 yra with the four spots 
and the six measured distances in different colours. 



Figure A2. Median flux measured in the wavelength calibration 
processed as science data as a function of the DRH wavelength 
(dashed black curve) and after correction (red plain curve). The 
dashed vertical lines show the wavelength of the laser lines that 
are used for the wavelength calibration of the IFS. 


For the wavelength calibration, the detector observes a flat 
source illuminated by four laser lines through the IFU, cre¬ 
ating four spots with precisely known wavelengths for each 
lenslet. In the DRH recipe, for each lenslet, the position of 
these four spots is measured and a third order polynomial is 
fitted to determine the wavelength of all pixels correspond¬ 
ing to the lenslet. For the step (2) of our wavelength calibra¬ 
tion refinement, we make use of the wavelength calibration 
raw frame, which we process in exactly the same manner 
as the science data. The resulting (x,y. A) data cube shows 
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flat images with flux variations as a function of wavelength, 
which correspond to the presence or not of laser lines in the 
original raw frame. The median flux in each of the spec¬ 
tral channels of this cube is plotted as a function of the 
DRH-determined wavelength in Fig. A2. The 4 laser lines 
are clearly visible, but the peaks appear shifted with respect 
to the wavelength of the lasers, indicating some error or off¬ 
set in the DRH wavelength calibration. To correct for this, 
we implemented a Levenberg-Marquardt least-squares min¬ 
imization routine in which we perform the following steps: 

• We make a Gaussian fit of the position of the laser peaks 
as a function of the spectral channel index. 

• We attribute a wavelength to one of the channels (ref¬ 
erence channel) and we extrapolate the wavelength of all 
the other channels using the scaling factor discussed previ¬ 
ously. For the reference channel we take as first guess the 
wavelength determined from the DRH. 

• We determine the expected wavelengths at the position 
of the peaks determined in the first step by interpolating 
(linearly) the newly defined wavelength calibration. 

• Finally, we measure the difference between these ex¬ 
pected wavelengths and the known wavelengths of the lasers. 

The only free parameter is the wavelength of the reference 
channel, and the minimization is performed on the maxi¬ 
mum of the difference (in absolute value) calculated in the 
last step. The median flux as a function of the corrected 
wavelength is also plotted for comparison in Fig. A2. This 
time the peaks properly match the expected wavelengths, 
and the difference is between 0.37 and 0.85 nm, while it was 
between 1.7 and 19.4 nm with the DRH wavelength calibra¬ 
tion. Overall, combining the errors on the determination of 
the scaling factor and the wavelength of the reference chan¬ 
nel, we estimate an error on the wavelength calibration of 
the order of 1.4 nm at most. 

Our procedure assumes that the position of the spectra 
on the IFS detector, and hence the wavelength calibration, 
is stable between the time when the science data is taken 
and the time when the wavelength is calibrated. SPHERE 
is located on one of the Nasmyth platforms of VLT Unit 
Telescope 3, which is mechanically very stable. The main 
temporal evolutions of the calibrations that are observed in 
SPHERE are due to bench temperature variations. For the 
IFS, the position of the spectra evolves by -^0.3 pixel/deg 
(internal ESO monitoring). During the Sirius observations, 
the temperature of the IFS was stable at 15.9°C, and for 
the wavelength calibration taken in the morning the tem¬ 
perature was 15.6°C, resulting in a shift of 0.09 pixel of the 
spectra. We consider this variation to be negligible and do 
not attempt additional corrections for this effect. 

A3 Scaling of the off-axis proxy PSFs 

As explained in Section 3.3, the off-axis PSF was slightly 
saturated in the first five spectral channels of the IFS, and 
the core of the PSF in the K12 Alters was saturated in IRDIS 
over a diameter of ^X/D. 

For the scaling the /3 Pictoris proxy PSF, we base our 
analysis on the first Airy ring of both PSFs, which is well 
within the linearity range for Sirius with ^3 000 counts on 
average. To determine the scaling ratio, we use a fitting rou¬ 
tine with three free parameters: the shift in x, the shift in 
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Figure A3. Scaling factor necessary for the 0 Pictoris proxy PSF 
to match the photometry of Sirius A in IRDIS the K12 filter pair. 
The factor has been determined by fitting the proxy PSF on in¬ 
dividual DITs of the Sirius PSF (black diamond and associated 
error bars). The dotted line corresponds to the average of all val¬ 
ues, and the grey area to the average plus or minus one standard 
deviation. 

y, and the flux scaling between the Sirius and proxy PSFs. 
During the fitting, the proxy PSF is shifted, scaled, and sub¬ 
tracted to the Sirius PSF. For IRDIS, we evaluate the stan¬ 
dard deviation of the residuals in an annulus with an inner 
radius starting at 1.5 X/D, and an outer radius Rout that we 
vary between 3 and 6 X/D by steps of 0.5 X/D. Figure A3 
illustrates the average value and standard deviation of the 
scaling factor determined with the different values of Rout 
for IRDIS in the K1 filter. We perform this analysis on the 
three off-axis PSFs that were acquired for Sirius (30 NDIT 
for each). We adopt as the final scaling factor the average 
of all values determined for all values of Rout, 87.4 ±4.9. As 
a sanity check, we note that the difference of K band mag¬ 
nitude between jS Pictoris and Sirius is equal to 4.83, corre¬ 
sponding to a factor of 85.5 that is perfectly consistent with 
our determination. The uncertainty on the determination of 
the scaling factor translates to a photometric uncertainty of 
0.06 mag on the final PSF. 

For the IFS, a similar analysis is done for the five first 
and four last spectral channels in which the PSF is just 
slightly saturated or out of the linear range. In these chan¬ 
nels, the SNR of the Sirius PSF outside of the first Airy ring 
is not enough to extend Rout to more than 3 X/D so we fit 
the proxy PSF only with this value. The estimated photo¬ 
metric uncertainty is of ~0.2 mag for these channels. As an 
additional check, we compared the flux of the corrected off- 
axis PSF with the flux of the satellite spots in the star centre 
images, and they agree within the 0.2 mag uncertainty. The 
other channels are not saturated so no particular correction 
is applied to the PSF in these channels. 

A4 Centring of the IFS coronagraphic data 

In Section 3.4 we described how the stopped near-infrared 
ADCs affect the centring of the star behind the corona- 


MNRAS 000, 1-16 (2015) 















High-contrast imaging of Sirius A with SPHERE 15 



Figure A4. Spatially rescaled images in three spectral channels of the IFS taken from the same data cube. The red circles delimit the 
annulus where the residuals of the difference between the reference frame and the others are minimized in the centring procedure (see 
text for details). The green dashed circle shows the position and size of the coronagraphic mask. Each frame is centred with respect to 
the star. Because of the ADC tracking problem, the star centre and the coronagraphic mask centre do not correspond at all wavelengths. 
However, the speckle pattern and diffraction spikes of the telescope spiders match very well in all the channels. The optimization annulus 
is chosen to avoid the bright speckles at the edge of the coronagraph that vary a lot with wavelength because of the sub-optimal star 
centring, but also the aliasing at the border of the AO control radius. 


graphic mask, especially towards the shorter wavelengths 
in the IFS data. Fortunately, the DTTS ensures that the 
centring is accurate close to a wavelength of 1.5 ^m, where 
it measures the differential tip-tilt between the visible and 
the near-infrared. We verified that it is the case in our data 
using the star centre frames (see Section 3.4). We adopt 
as the reference wavelength the spectral channel number 35 
(A 35 = 1.58 /rm), which means that we consider that the 
star is well centred on the coronagraph in this channel, and 
that it remains stable throughout the observation, despite 
the ADC problem. We use as reference centre in this chan¬ 
nel the average of the three centres obtained from the three 
star centre frames. 

The second step of the analysis, for each IFS data cube, 
consists in using this reference channel to recentre all of 
the others. First, all channels are rescaled spatially using 
the wavelength-dependent scaling factor determined for the 
wavelength calibration in Appendix A2. The value of this 
scaling factor does not vary with time and can be deter¬ 
mined once and for all for a given observation using the star 
centre frames. For an accurate spatial rescaling, we used an 
FFT-based procedure, detailed in Appendix B, which per¬ 
forms zero-padding both in the direct and Fourier spaces. 
After the rescaling, the speckle pattern has the same size in 
all frames, and we can find the shift between any spectral 
channel and the reference channel. In practice, the optimal 
shift is determined by minimizing the residuals in the differ¬ 
ence between the reference channel and the other channels. 
The residuals are calculated in an annulus between 30 and 50 
pixels. This area has been chosen because it is within the AO 
control radius in all channels but avoids the varying effect 
of aliasing at the edge of the control radius, and it contains 
bright speckles that can be easily matched between channels 
while avoiding the bright residuals close to the coronagraph 
edge. These residuals are particularly strong and varying be¬ 
cause of the differential refraction that induces a decentring 
of the star behind the coronagraph with wavelength. Fig. A4 


shows three spectral channels from the same cube with the 
optimization annulus, and the position and size of the coro¬ 
nagraphic mask. The frames are centred with respect to the 
star, resulting in a visible motion of the coronagraphic mask 
position as a function wavelength. However, the pattern of 
instrumental speckles is stable and most speckles within the 
optimization annulus can be seen in all three channels. 


APPENDIX B: SPATIAL SCALING 
ALGORITHM DETAILS 

Accurate spatial rescaling of discrete images is very impor¬ 
tant for high-contrast imaging where techniques such as SDI 
(Racine et al. 1999) are heavily used for speckle noise at¬ 
tenuation. We detail below an FFT-based rescaling tech¬ 
nique that makes use of zero-padding in both the direct and 
Fourier spaces for accurate rescaling of discrete data. 

For conciseness we consider one-dimensional signals, 
and the results generalize readily to two dimensions. Given 
N samples of a signal i sampled with a pitch pa , the quantity 
i{k.pi,) given by: 

JV-l 

i{k.p„) = ^ (Bl) 

1=0 

is the discrete Fourier transform of i at frequency k.p,^, where 
Pi, is the pitch in Fourier space. Because of the use of the 
fast Fourier transform (FFT), the pitch in Fourier and direct 
spaces are related by 

PaPv = 1/N, (B2) 

which means that the support Npa is the inverse of the 
Fourier pitch. 

If we modify the number of pixels in direct space by 
adding on each side of the table some null pixels before per¬ 
forming the FFT, equation B2 becomes PaPv' = l/fV^ where 
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N' is the new signal size after zero-padding in direct space, 
and p^i the new pitch in Fourier space. 

If we proceed similarly in Fourier space by adding null 
pixels after performing the FFT, we finally obtain for the 
new spatial sampling in direct space after an inverse FFT: 


P^’> = {N'/N”)p^, (B3) 

where N" is the new signal size after zero-padding in 
Fourier space. We have thus zoomed the image by a fac¬ 
tor 2 = {N”/N'). Note that if we only zero-pad the image 
in Fourier space, as it is often done, then N' = N, which 
severely constrains the zooming factor to 2 = N” /N. 

Considering a coronagraphic image Ai that we want to 
resample at a smaller wavelength Aq , we need to satisfy the 
following equation: 


N' _ Al 
^ Ao ■ 


(B4) 


Having two degrees of freedom N' and N” allows us, 
in practice, to satisfy this equality with a good precision 
(typically 10“^) for reasonable values of N' and N” for all 
wavelengths of the IRDIS or IFS instruments. 


This paper has been typeset from a TFX/IATI^ file prepared by 
the author. 
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